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ABSTRACT 

We demonstrate that measurable vertical structure can be excited in the accretion 
disc of a close binary system by a dipolar magnetic field centred on the secondary 
star. We present the first high resolution hydrodynamic simulations to show the ini- 
tial development of a uniform warp in a tidally truncated accretion disc. The warp 
precesses retrogradely with respect to the inertial frame. The amplitude depends on 
the phase of the warp with respect to the binary frame. A warped disc is the best 
available explanation for negative superhumps. 

Key vifords: 

accretion, accretion discs — instabilities — hydrodynamics — methods: numerical — 
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1 INTRODUCTION 

Much of the variability of cataclysmic variables can be 
traced to instabilities in the accretion disc. Thermal-viscous 
instability leads to dwarf nova outbursts, and tidal instabil- 
ity leads to superhumps and superoutbursts. More recently, 
people have been giving thought to warp-type instabilities. 
Pringle (1996) showed that radiation from the central accret- 
ing object will exert an asymmetric force on either surface 
of the disc that can be sufficent to warp it. However, rel- 
atively extreme conditions are required to generate a warp 
in this way (Ogilvie & Dubus 2001). In this paper we show 
numerically that stellar magnetic fields can warp the disc of 
a cataclysmic variable. In particular we shall focus on the 
effects of the secondary (donor) star's magnetic field. 

The light-curves of many cataclysmic variables exhibit 
periodicities close to but distinct from the orbital period. 
These signals are categorised as positive superhumps if they 
are longer than orbital and negative superhumps if they are 
less than orbital. Both are thought to emanate from the 
disc with the well established explanation for the former be- 
ing that they are the visible consequences of an eccentric 
disc precessing with respect to the binary's tidal field (see 
O'Donoghue 2000 for a review). Negative superhumps have 
not been studied in nearly as much detail. The observations 
of negative superhumps are summarised in Patterson (1999) 
in which 11 systems are listed. With periods typically a few 
per cent less than orbital, a key feature of negative super- 
humps is that they appear to be independent of positive 
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superhumps. They have been observed both in the absence 
of and in conjunction with their better studied brethren. 

Patterson et al. (1993) suggested that negative super- 
humps occur when the accretion disc is tilted to the binary 
plane. Just as the line of nodes of the Moon's orbit about 
the Earth regresses, so too the intersection of the disc and 
binary planes precesses retrogradely. The effect is present in 
the purely dynamical problem and is analogous to the eccen- 
tric disc case. A test particle in circular orbit around a single 
star has three natural frequencies, that of its azimuthal mo- 
tion (uj), and those of any radial or vertical motion. Around 
a single star these frequencies are identical and any pertur- 
bation to a circular orbit causes the particle to orbit in an 
ellipse. For a particle in orbit around a star in a binary we 
find that the tidal field breaks the symmetry of the three 
periods. The radial frequency becomes slightly lower than 
(ui) and the vertical frequency slightly higher. As a result, 
the axis of any perturbation to the orbit will precess. An- 
alytical calculations (e.g. Papaloizou & Terquem 1995) and 
numerical simulations (e.g. Larwood et al. 1996; Wood et al. 
2000) showed that an initially rigidly tilted disc in a binary 
system will precess retrogradely, approximately as a rigid 
body, as a result of the gravitational torque acting on the 
disc. 

The major stumbling block for the warped disc hypoth- 
esis has been the lack of a viable mechanism for generating 
it. Three possibilities can be explored. A warp might have 
developed as a result of an impulsive torque or out of the 
plane initial condition, and then precessed freely (see e.g. 
Larwood et al. 1996). A second possibility is that the warp 
developed from a co-planar disc via an instability. Lubow 
(1992) and Pringle (1996) discussed dynamical and radia- 
tive disc instabilities respectively. The third possibility, the 
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one explored in this paper, is that the warp is excited by 
continual forcing out of the plane. 

Both impulsively and resonantly generated waxps would 
be forced back into the plane as the disc evolved viscously. 
Murray & Armitagc (1998) presented a series of calcula- 
tions suggesting that the dynamical and radiative instabil- 
ities were too weak to overcome viscosity. Whilst continual 
forcing could overcome the problem of viscous flattening, 
the long term disc response will not be free precession but 
be with the same period as the forcing. 

The inspiration for this paper was the idea espoused by 
Terquem and Papaloizou (2000), that a warp might be in- 
duced in a circumstellar disc by the stellar magnetic field. 
Should the axis of the magnetic field be misaligned with 
the binary rotation axis, then the magnetic field will exert 
a torque in the plane of the accretion disc. Terquem and 
Papaloizou were interested in T Tauri stars, however the 
principle could equally bo applied to a cataclysmic variable 
with a moderately magnetised accreting white dwarf. Lai 
(1999) described in detail the mechanisms by which the stel- 
lar magnetic field would couple with the disc. Again however 
his attention was directed more at X-ray accreting pulsars 
and not at cataclysmic binaries. 

We carried out an initial series of smoothed particle 
hydrodynamics simulations and did indeed find that verti- 
cal structure was induced in the disc. However this struc- 
ture co-rotated with the accreting star and could not pos- 
sibly generate the periods characteristic of negative su- 
perhumps. We will describe these simulations in a sub- 
sequent paper (Wynn, Murray & Chakrabarty in prepa- 
ration). Pringlc (private communication) suggested we al- 
low the donor rather than the accretor to be magnetised. 
Whereas the primary star in a cataclysmic variable typically 
rotates several times faster than the binary, the secondary 
is tidally locked and hence co-rotates with the binary. This 
paper describes a series of simulations with which wo sought 
to determine whether the resultant accretion disc structure 
was approximately stationary in the binary or in the inertial 
frame. 

Several authors have already performed detailed analy- 
sis of the suitability of the smoothed particle hydrodynam- 
ics (SPH) technique for modelling three dimensional struc- 
ture in accretion discs. Murray (1997) and Yukawa, BofRn 
& Matsuda (1997) demonstrated that tidally excited spi- 
ral structure is resolved in well constructed two dimensional 
and three dimensional SPH simulations respectively. Lar- 
wood et al. (1996) modelled the precession of tilted discs 
in close binaries. They followed the long term evolution of 
discs tilted at 45 and 90 degrees to the binary plane us- 
ing 17500 particle simulations. Nelson & Papaloizou (1999) 
studied the propagation of bending waves in accretion discs. 
They specifically tested SPH simulations against the linear 
theory of warped discs and found good agreement. Calcula- 
tion resolution ranged from 20000 to 102000 particles. They 
later (Nelson & Papaloizou 2000) studied the warp produced 
by the Bardeen-Pcttcrsou effect using simulations ranging 
up to 200000 particles in size. Wood, Montgomery & Simp- 
son (2000) used simulations with 25000 particle resolution 
to follow the precession of a prc-tilted disc in a close binary. 

The calculations described later in this paper involve 
over 500000 particles. At this resolution we are not limited 
to resolving the gross features of the warp itself. We can 



also determine the the vertical structure caused by the tidal 
forces of the secondary, and then follow the interaction be- 
tween the two. 



2 MAGNETISED SECONDARIES 

The secondary stars in cataclysmic variables have rotational 
periods which are tidally locked to the binary's orbital pe- 
riod, and are typically of the order of a few hours. The 
connection between rapid stellar rotation rates and chro- 
mospheric, and by implication magnetic, activity has long 
been established (eg. Kraft, 1967; Skumanich, 1972). This 
connection is often invoked as an argument in favour of the 
magnetic nature of the rapidly rotating secondary stars in 
CVs. Direct observational evidence for magnetic activity re- 
mains elusive however. Some of the stongest evidence comes 
from observations of star spots (and their close association 
with magnetic activity) in the CV ST Leonis Minoris (How- 
ell et al., 2000) and the pre-CV V471 Tau (Applegate & 
Patterson, 1987) as well as in a number of other types of 
close binary systems. The evidence for solar-type activity 
cycles in the secondary stars of CVs is discussed by Bian- 
chini (1990). Any such cycle would imply the presence of 
significant magnetic fields. The existance of these magnetic 
fields has strong theoretical support. Angular momentum 
loss by magnetic braking of the secondary star is the gen- 
erally accepted mechanism driving the orbital evolution of 
CVs above the orbital period gap (Verbunt & Zwann, 1981; 
Mestel & Spruit, 1987). Moreover, theories which attempt 
to explain the sychronism between the binary orbital period 
and the rotational periods of magnetic white dwarfs in CVs 
require secondary magnetic moments of the order /i2 ~ 10^^ 
- 10^"* G cm=^ (eg. King, Whitehurst & Prank, 1990). 

For the purposes of the simulations presented in this pa- 
per we assume that the secondary stars in CVs have a signif- 
icant dipolar magnetic field, and that the magnetic moment 
is similar to that quoted above. This implies that the (un- 
perturbed) field strengths close to the edge of the accretion 
disc are of order ~ 10^ G. As the gas in the accretion disc 
orbits the white dwarf in the binary it will twist the field 
of the secondary star and induce currents within the disc. 
The interaction of these currents and the imposed magnetic 
field gives rise to a variety of forces. King (1993) developed 
a drag prescription to parametrise the effects of these mag- 
netic forces on accreting gas in the AM Herculis sub-class of 
cataclysmic variables. Pearson, Wynn and King (1997) later 
made use of this prescription to study the effects of the az- 
imuthal force arising from the interaction of radial surface 
currents and the vertical component of the magnetic field on 
the structure of a two dimensional disc in the orbital plane 
of the binary. This study showed that disc-field interaction 
caused the disc to lose angular momentum (to the secondary 
star via the field lines) causing the disc radius to decrease, 
and enhancing the mass accretion rate. Moreover, resonance 
phenomena similar to those in the SU UMa systems caused 
the disc to become eccentric. In this paper we extend the 
analysis of the effect of the magentic field on the accretion 
disc to three dimensions. 

Lai (1999) made a detailed analytic study of the torques 
to which a disc is subjected to when an external stellar 
field is imposed upon it. In the general case he found that 
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a warping torque resulted from the interaction of the az- 
imuthal component of the stellar magnetic field and the ra- 
dial surface current on the disc induced by the twisting of the 
threaded vertical field component. However, the details of 
the disc-field interaction remain unclear, in particular there 
is little agreement on the range and strength of the magnetic 
stresses over the disc. In order to examine the response of 
a fiuid disc to an applied warping torque we implemented a 
simplified version of the magnetic interaction in which the 
gas within the disc is subject to an acceleration of the gen- 
eral form 



(1) 



where v and Vf are the velocities of the material and field 
lines respectively, and the suffix _L refers to the velocity 
components perpendicular to the field lines. The parameter 
k - 1/t mag reflects the details of the plasma-magnetic field 
interaction, and for the purposes of this paper we adopt the 
simple scaling k oc B^{r, 9, 0) and assume that the magnetic 
field has a purely dipolar form. We can estimate the magni- 
tude of k in the vicinity of the accretion disc in a cataclysmic 
variable with a magnetic secondary star from the compari- 
son of (|l|) and the calculations of Lai (1999). Equating the 
force per unit area in each case we have 



2 



(2) 



where r2 is the distance from the secondary star and E is the 
surface density of the accretion disc. Typically r2 ^ 10^° cm, 

- 100 km 



10^* G cm^ E ~ 100 g cm"^ and 



M2 

s~^ in a CV giving ^(10^° cm) 



10" 



[V-V/ 

„-i 



This simplistic 



form of the magnetic force preserves all of the components 
required to induce a warp within the disc and has a mag- 
nitude of the order of that expected of magnetic stresses 
induced on the accretion disc by a magnetic secondary in a 
CV. 



3 CODE ASPECTS 

We use a smoothed particle hydrodynamics code that has 
been described extensively elsewhere (Truss et al. 2000, Mur- 
ray 1996). The code uses operator splitting to separate the 
calculation of the gravitational and magnetic forces from 
that of the gas-dynamical processes. The code has been par- 
allelised using OpenMP to run on shared memory paral- 
lel supercomputers. The viscosity is slightly modified from 
Murray (1996). We are free to choose the length scale in the 
viscosity term so we set it to equal the disc scale height, 
and so are able to reproduce a Shakura-Sunyaev a viscosity. 
The viscosity term also produces a bulk viscosity of similar 
magnitude in regions where the velocity field is significantly 
divergent. A variable smoothing length h is used. To en- 
sure the disc is everywhere resolved h must be less than 
the scale height whilst maintaining an adequate number of 
neighbours. 

We use an isothermal equation of state. In other words 
we assume any viscously dissipated energy is instantly radi- 
ated away. 

The magnetic interaction is iinplemented by including 
the acceleration given by equation □ along with the gravita- 
tional terms in the operator split equations of motion. 





Figure 1. Surface density map of the initial state. Note the finely 
resolved tidal spiral arms. There are 536033 particles in the disc. 
Midway into the disc there is a ring where the surface density is 
a maximum. This is an artefact from the disc's creation. The key 
to the right shows the colour coding from low density (red) to 
high density (black). 
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Figure 2. Silhouette of the planar (unwarped) disc as seen by 
an observer in the binary plane. Azimuth —it is along the binary 
axis between the stars, and then increases anticlockwise round 
the disc. 



As in previous papers we have normalised units so that 
the binary separation d, the total system mass M, and the 
binary angular velocity Qb are all equal to one. 
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Figure 4. Vertical structure of tlie disc at times t=5,6,7,8,9 and 10 (from top left to bottom right). Blue liues indicate structure below 
the disc plane, green shows regions that are symmetric about z = whilst yellow and red show regions that are progressively more lifted 
above the disc plane. The panels in this figure correspond to those in figure H. 
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Table 1. Simulation parameters. The magnetic dipole on the 
secondary is oriented so as to be (j>d to the vertical and 0^ prograde 
of the binary axis. 



parameter 


symbol 


value 


mass ratio 


<? 


0.4 


sound speed 


c 


0.05 d^i, 


viscosity parameter 


a 


1.0 


magnetic parameter 


k 


10-5 (10-2 s-i) 


dipole declination 


4>d 


20 deg 


dipole azimuth 


9d 


15 deg 



4 CALCULATIONS 

The system and code parameters for the calculations are 
listed in Table ^. We constructed a three dimensional initial 
state for the calculation by running a two dimensional non- 
magnetic simulation to equilibrium. The particles were then 
replicated in z to produce a three dimensional disc which was 
then followed for several orbits to allow vertical oscillations 
to die out. The resulting disc comprised 536033 particles. It 
is shown in plan view, shaded according to surface density, 
in figure 0. Tidally inspired spiral shocks are clearly visi- 
ble. We see no evidence that spiral shocks are suppressed 
or any weaker in three dimensions than in two dimensional 
calculations. 

We also show the silhouette of the initial disc (figure ^ 
as would be seen by an observer on the disc plane. This 
figure was obtained by plotting all the particles on the 9 — z 
plane. The two spiral arms produced 50 per cent increases 
in the disc height which should be observable in X-ray and 
EUV observations of eclipsing systems (see e.g. Billington 
et al., 1996). 

As we shall see below, the starting disc was not in equi- 
librium with the mass addition stream from the secondary 
star. However the time-scale for mass evolution is much 
longer than that for magnetically induced evolution and this 
did not affect our results. 

We completed two simulations; a reference calculation 
with no magnetic field, and a second in which the magnetic 
field was added instantaneously. The dipole orientation 

jl = sin((?id)(cos(6lo + ^bt) i + sin(6lo 4- ^bt)i) + cos((;/>d) k(3) 

where the angular speed of the binary fii, = 1. i,j and k are 
the cartesian unit vectors. In our coordinates the disc lies 
in the x-y plane with the angular momentum vector of the 
binary in the positive z direction. 

Mass transfer from the secondary was incorporated into 
the simulations by adding single particles from L\ at inter- 
vals A t = 0.01 Porb- As in previous simulations (e.g. Murray, 
1996), the initial velocity of injected particles was taken from 
Lubow & Shu (1977). The disc was not in equilibrium with 
the mass transfer stream which was only incorporated into 
the simulation as a useful tracer of the secondary magnetic 
field. Gas fiow on the secondary star was not resolved in the 
calculations. 

The calculations described here (and their antecedents) 
took the order of one week of wall time on 8 processors of 
an SGI Origin 2000 or 3000. 
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Figure 5. Strength of the warp as a function of time when the 
secondary is magnetised (solid line) and unmagnetised (dotted 
line). Both curves have been normalised by multiplying by 10'^. 
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Figure 6. Average disc radius for the magnetised (solid line) and 
unmagnetised (dotted line) simulations. 



5 RESULTS 

In this section we show that the drag on an accretion disc, 
exerted by a tilted magnetic dipole anchored to the sec- 
ondary star, results in a warp. Figure ^ is a sequence of six 
disc silhouettes, each generated in the same manner as fig- 
ure ^ The snapshots show the disc at equally spaced times 
from t = 5 — 10 (in our units Porb = Stt). Whilst the disc 
is approximately symmetric about its midplane in the first 
three panels, by the fourth {t = 8) it is warped. 

We have plotted six plan views of the disc vertical struc- 
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ture in figure corresponding to the snapshots in figure ^ 
Each image is constructed by interpolating particle z (aver- 
aged over all particles in a region) to a grid. Blue regions of 
the disc are predominantly below the z = plane. Green re- 
gions are approximately symmetric about the midplane, and 
yellow-red regions are more extended above z — 0. The mag- 
netic torque deflects the stream of gas from the secondary 
above the plane (redwards) . Contours have also been drawn 
to aid the eye. Figures ^ and ^ therefore show the disc's 
initial response to an instantaneously applied out of plane 
force. 

Less than one binary period after the magnetic force 
was applied, the disc warp was a coherent global structure 
(panel 1). However one sixth of a binary period later (panel 
2) the warp was no longer apparent. It reappeared again, 
much stronger, in panels three and four. In each image the 
warp is a disc- wide structure, with the spiral density waves 
in the initial disc having been rapidly diminished by the 
warp. 

If one follows, for example, the yellow-red zone from 
panel to panel one sees that the warp has precessed retro- 
gradely in the binary frame. The disc precession with respect 
to the inertial frame occurred on a much longer time scale 
so quantitative analysis of the calculations was required. 

We Fourier transformed the disc's vertical structure in 
both azimuth and time, as in Murray & Armitage (1998). 

The sine transform of the vertical displacement, 

p=i 

where Op is the angular position of particle p (measured 
in the inertial frame) and A*' is the total number of parti- 
cles. The cosine transform is defined similarly and the mode 
strength is set equal to the root mean square of the two. The 
mode strength has units of length, and so a large response 
would be of the order of the disc height (0.01 d). As a check, 
we returned to the initial disc (as shown in figure |l|) and 
applied a uniform tilt of five degrees to it. We then evolved 
the tilted disc for one orbital period, and applied equation ^ 
We found ~ 0.02 which agreed with the average vertical 
displacement of the particles. 

Figure ^ overlays the warp mode strength for the two 
calculations. For convenience only we have normalised both 
curves, by a factor 10^. No warp developed in the refer- 
ence calculation, and the corresponding curve (dashed line) 
is noisy but with typical magnitude of only 10~* or a factor 
two hundred smaller than for the five degree warp. 

When the magnetic field was switched on in the second 
calculation, the disc immediately responded (solid line). In- 
stead of increasing to a steady value however, the warp mode 
strength oscillated with a period P„ — 0.987 ± 0.003 Porb as 
the warp slowly precessed retrogradely in the inertial frame, 
and the structure of the warp changed as it moved relative 
to the binary. High time resolution animations of the disc 
motion were also made that confirmed this to be the pre- 
cession period. Typical amplitudes of 5*^ for the magnetic 
calculation were of the order of 10"'^, about a factor ten or 
twenty less than for the tilted disc. These amplitudes then 
agree with the vertical structure shown in figure ^. 

There was a secular decline in warp amplitude. Now 
in these simulations we had a = 1 and disc scale height 
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Figure 7. Disc mass for the magnetised (solid line) and unmag- 
netised (dotted line) simulations. 
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Figure 8. Disc dissipation for the magnetised (solid line) and 
unmagnetised (dotted line) simulations. 

H ~ 0.04r. One can estimate the time scale for viscous forces 
to fiatten a warp tf ~ a{r/H)^ . At the disc outer edge 
tf ~ 300, not entirely inconsistent with the decay time scale 
apparent in figure 0, particularly when one considers the 
disc shrinkage over the course of the magnetic calculation. 
Reducing the effective a of the disc would result in more 
rapid damping of the warp. 

The applied magnetic torque swiftly removed angular 
momentum from the outer disc, causing it to shrink. Figure ^ 
shows the evolution of the average disc radius (Rd), which we 
define to be the radius of a circle with the same area as the 
disc. The figure clearly shows that the disc shrinks after the 
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application of the magnetic field. This was predicted by the 
earlier treatment of Pearson, Wynn and King (1997) and is 
a consequence of the relation tvisc(-Rd) ~ imag(-Rd)- Figure ^ 
shows the mass declined more rapidly for the torqued disc 
than for the reference calculation. Energy dissipation (which 
can be equated to bolometric luminosity for the disc) due to 
viscosity in the two discs is plotted in figure |^. The sharp rise 
in dissipation correlates with an accelerated rate of accretion 
onto the white dwarf. It is no surprise therefore that the 
warp strength declined in figure ^. 

Figures ^ and ^ show the changes in silhouette and 
plan view respectively over the time period t = 22 — 27. The 
retrograde precession of the warp can clearly be seen in fig- 
ure hd. The variations in warp amplitude are also apparent. 
The warp almost disappeared when it was 180° out of phase 
with the magnetic torque (in between panels 3 and 4). Even 
then however the disc was not flat, but the vertical structure 
that was present at that time was symmetric about the mid- 
plane (figure ^ panels 3 and 4) . Half a binary period later 
(between panels 6 and 1 as the structure changed cyclically) 
the warp was reinforced by the magnetic torque and the disc 
was most distorted. In panel 1 for example the scale height 
varied by as much as 30 % around the circumference of the 
disc. Such variations are sufficiently large to be observable 
in nearly edge on systems. Again we refer to the Billington 
et al. (1996) analysis of ultraviolet superdips in OY Car. 

A simple toy model, constructed by the anonymous ref- 
eree, is very helpful in understanding the results shown in 
figures §,|and0. Consider a single annulus of the disc, and 
neglect pressure and viscosity. If the inclinations are small, 
then the tilt of the annulus can be described by a variable 
W(/:) = l^-\-ily (where x and y are spatial coordinates mea- 
sured in the binary frame). W(t) obeys a linear equation of 
the form 
dW 

— =iLjpW-7(W- Wb), (5) 

where Wp < is the tidal precession rate in the binary frame, 
and 7 represents the effect of the magnetic drag. In the ab- 
sence of precession, the drag forces would bring the disc tilt 
W to coincide with Wb (determined by the orientation of 
the secondary magnetic field) on a timescale I/7. Starting 
with zero initial tilt, we have 

W(f)=WB [1 -exp(iLjp -7)t]. (6) 

7 - 

The oscillations in tilt amplitude in figure ^ are neatly ex- 
plained as |W| is oscillatory with frequency ujp. Equation ^ 
also shows that the warp's precession is a transient response 
to the initial kick of the magnetic field. The steady state 
disc response will be a disc warp that is stationary in the 
binary frame. 

Negative superhumps have not been subject to the same 
degree of scrutiny as their positive counterparts. We are con- 
fident that positive superhumps arise in a precessing disc 
(see most recently Rolfe, Haswell & Patterson, 2001). Obser- 
vations of negative superhumps are limited to photometry, 
with no information regarding either temperature or loca- 
tion of the emitting region. Patterson et al. (2001) describes 
a persistent negative superhump observed in the nova-like 
variable V751 Cygni, and states that "... many nova-like 
variables in their bright state and with Porb = 3 — 4 hr show 
negative superhumps" . They do not show the close associa- 



tion with superoutbursts that normal superhumps do, and 
they do not appear in very low mass ratio systems (i.e. those 
below the period gap). 

Positive superhumps are caused by a tidal resonance 
which can only act upon an accretion disc if the ratio of 
secondary to primary stellar masses q ^ 1/3 (Lubow 1991; 
Murray, Warner & Wickramasinghe 2000). Lubow (1992) 
showed that a second resonance that acts to tilt the disc, 
almost exactly overlays the eccentric resonance. The tilt res- 
onance is however much weaker, and simulations described 
in Murray & Armitage (1998) failed to produce a disc warp, 
but with 30 000 particles in the calculations any detection 
of a warp would have been marginal. 

It is worth revisiting the problem, making use of the 
much higher resolution now available to us (10® particles) 
and also considering a higher mass ratio. The strength of the 
tilt resonance goes as and mass ratios of observed neg- 
ative superhumpers are closer to 1/3 rather than the 3/17 
used in Murray & Armitage. The calculations presented in 
this paper give us confidence that smoothed particle hydro- 
dynamics simulations are an appropriate tool for studying 
the evolution of warp modes in discs, however the warps are 
generated. 



6 CONCLUSIONS 

We have followed the evolution of a warp in a close binary ac- 
cretion disc. An initially planar disc was subjected to a force 
that modelled the effects of an inclined magnetic dipole cen- 
tred on the mass donor star. The magnetic field was assumed 
to be fixed in the binary frame. 

Other authors have modelled the evolution of discs with 
pre-existing warps, or of discs that were misaligned with the 
binary plane. However these are the first hydroydnamic sim- 
ulations to follow the initial development of a warp, and the 
first to show how the warp structure varies as it precesses. 

The disc's initial response to the inclined field was com- 
plicated by spiral shocks excited by the secondary's tidal 
field. However the spiral shocks are strongly damped by the 
developing warp. In turn the warp structure rapidly simpli- 
fied and steady, uniform retrograde precession (in the iner- 
tial frame) occurred. 

The amplitude and structure of the warp depended 
upon its phase with respect to the perturbing force (and 
hence of the binary). The warp was maximum when it 
was reinforced by the magnetic force. However 180 deg 
later, when the warp and the force were anti-phased, the 
disc was almost flat. We measured the period of the mo- 
tion of the warp with respect to the binary to be Pw = 
0.987 ± 0.003 

Porb ■ 

We suggest therefore that negative superhumps could 
be explained not simply via the retrograde precession of a 
warped disc, but in terms of changes in the discs vertical 
structure that occur on the same period, Pw. Because of the 
very limited observations of negative superhumps we can but 
hypothesise as to the precise location of the luminosity vari- 
ations. One possibility is the changing disc aspect presented 
to the accretion stream as the warp precesses. Alternatively, 
negative superhumps may be the result of radiation from 
the central source being reprocessed over a changing disc 
surface. 
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Figure 9. Disc silhouette at times t=22,23,24,25,26 and 27. Panels 4 ans 5 show an almost flat disc and correspond to the minimum in 
warp strength (see figure ft). 
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In subsequent papers already in preparation we shall 
discuss both the evolution to steady state of discs in mag- 
netised systems, and the influence of magnetised accretors 
on disc structure. 
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